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Abstract 

The purpose of this article is to present the thermodynamics of the pseudospin 
electron (PE) model in the case of the different type interactions between 
pseudospins. First, we provide an overview of the results of works which deal 
with the theoretical investigation of the PE model with the inclusion of the 
direct pseudospin-pseudospin interaction (but without the electron transfer). 
Second, we present the results of the investigation of the model in the case 
of the absence of the direct pseudospin-pseudospin interaction and Hubbard 
correlation, when interaction between pseudospins via conducting electron is 
done. 
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Introduction 

Pseudospin-electron (PE) model is one of theoretical models which con- 
siders the interaction of electrons with local lattice vibrations where an 
anharmonic variables are represented by pseudospins. The theoretical 
investigation of the PE model is an enduring subject of interest at the 
quantum statistics department. 

The model is used to describe the strongly correlated electrons of 
CuC>2 sheets coupled with the vibrational states of apex oxygen ions 
Oiv (which move in the double-well potential) in YBaCuO type high- 
T c superconductors (HTSC) [ B]. Recently a similar model has been 
applied for investigation of the proton-electron interaction in molecular 
and crystalline systems with hydrogen bonds [ ^] . 

The model Hamiltonian is the following: 



H = E H * + E *« ^ - \ E J v 3 ' s !> w 

i ijcr ij 

Hi = U n^nn - A*E nia + g E n i° 3 l ~ ^Sf , 



where the strong single-site electron correlation U, interaction with the 
anharmonic mode (g-term), and the energy of the anharmonic potential 
asymmetry (/i-term) are included in the single-site part. Hamiltonian ([!]) 
also contains terms, which describe electron transfer Uj and direct inter- 
action between pseudospins Jy . The energy of the electron states at the 
lattice site is accounted from the level of chemical potential /i. 

In the case of strong coupling (g 3> t) and correlation (U ^> t) 
the perturbation theory can not be applied. It is reasonable to include 
these one in zero order Hamiltonian (single-site Hamiltonian Hi). Its 
eigenfunctions are build of the vectors \n^, n^, Sf), which form the full 
unit cell state basis [ [| : 

|1> = |Q,0, |>, |2) = |l,l,i) ) |3> = |0, 1, |>, [4> = (1,0, |>, (2) 

|1> = \0A~), |2> = |1, 1, — 1>, 13} = [0, 1, ~|>, |4> = |1,0,-|>- 

In the early works the main attention at the investigation of this 
model has been paid to the examination of electron states, effective el- 
ectron-electron interaction, to the elucidation of additional possibilities 
of occurrence of superconducting pair correlations. On the basis of PE 
model a possible connection between the superconductivity and lattice 
instability of the ferroelectric type in HTSC has been discussed [ El O. 
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A series of works has been carried out in which the pseudospin (S Z S Z ), 
mixed (S z n) and charge (nn) correlation functions were calculated. It 
has been shown with the use of the generalized random phase approxima- 
tion (GRPA) in the limit U — > oo [ ||, 0], that there exists a possibility of 
divergences of these functions at some values of temperature. This effect 
was interpreted as a manifestation of dielectric instability or ferroelec- 
tric type anomaly. The tendency to the spatially modulated charge and 
pseudospin ordering at the certain model parameter values was found 
out. The analysis of ferroelectric type instabilities in the two-sublattice 
model of high temperature superconducting systems has been made [ [5J . 
The influence of oxygen nonstoichiometry on localization of apex oxygen 
in YBa2Cu307_2; type crystals is studied in the work [ ||. 

The purpose of this article (lecture) is to present the thermodynamics 
of the PE model in the case of the different type interactions between 
pseudospins. First, we provide an overview of the results of works which 
deal with the theoretical investigation of the PE model with the inclusion 
of the direct pseudospin-pseudospin interaction (but without the electron 
transfer (tij =0)). Second, we present the results of the investigation of 
the model in the case of the absence of the direct pseudospin-pseudospin 
interaction and Hubbard correlation (Jy = 0, U = 0), when interaction 
between pseudospins via conducting electron is done. 



1. Direct interaction between pseudospins. 
1.1. Ferroelectric type interaction. 

The work [ || is devoted to the study of the PE model in the case of zero 
electron transfer (ty = 0). The direct interaction between pseudospins 
is taken into account. It is supposed to be a long-ranged that allows one 
to use the mean field approximation (MFA). In this approximation, the 
model Hamiltonian has the following form: 

H= Hi+ — Jn 2 , H l =-n ^ n^+Un^nn +g ^ n ia Sf - (h+Jn) S? . 

i a a 

The interaction (Jy ~ J/N) is taken as the ferroelectric type one; 
the order parameter rj = (S z ) does not depend on the unit cell index. 

Grand canonical potential and partition function of the model, cal- 
culated per one lattice site are equal to 

^- = -T\nZ l+ l -Jn 2 , (3) 



Zi = 2 



cosh^/i„+e- /3(c/ - 2 ' j) cosh/3(fc n -s)+2e ft * cosh/3 (hn~] 
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where h n = h/2 + Jr]/2. 

Then, all thermodynamic properties can be derived from the pre- 
sented formulae (||). 

The average number of electrons is determined as follows [ ^ : 

The equation for the order parameter is obtained from the thermo- 
dynamical relation [ ||] : 

For the investigation of equilibrium conditions one should separate two 
regimes: fi = const and n — const. We would like to note that hereafter 
we shall demonstrate the results of numerical investigation which show 
the main features of the considered model. 



1.1.1. /i = const regime. 

The ll = const regime corresponds to the case when, for example, charge 
redistribution between the conducting sheets CuC>2 and other structural 
elements is allowed. For this regime the equilibrium is defined by the 
minimum of the grand canonical potential that form an equation for 
pseudospin mean value (^) and expression for (n) (^J). 

At some regions of h values the held dependencies n(h) possesses S - 
like behaviour Fig. [j] (the first order phase transition with the jump of 
order parameter and electron concentration take place at the change of 
field h). The phase transition point is presented by a crossing point on 
the dependence fl(h). 

The phase diagram T c — h is shown in Fig. ||. One can see that with 
respect to Ising model the phase coexistence curve is shifted in field and 
distorted from the vertical line and hence the possibility of the first order 
phase transition with the temperature change exists in PE model. 

The phase diagram fi — h (Fig. |^) indicate stability regions for states 
with r) = ±1/2 (U < g and T = 0). The form of diagram depends on 
the relation between U and g parameter values [ |^] . Transitions between 
regions r) <-> |p), f) <^ \p) lead to the change of the average number of 
electrons only. At transitions \r) <-> \f) the flipping of pseudospin takes 
place, and at \r) <-> \p) (r ^ p) both processes occur. 
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Hence, the possibility of the first order phase transition with the 
change of field h and/or chemical potential /i takes place and is shown 
by thick line on phase diagram in Fig. |J 



1.1.2. 



const regime. 



0<n<l 



1/2-n 













In the regime of the fixed electron concen- 
tration value the equilibrium is determined 
by the minimum of free energy F = £l+/j,N 
and form a set of equations (||) , (||) for the 
chemical potential and order parameter. 

Typical example of the rj(h) depen- 
dence is shown in Fig. 3 which corre- 
sponds to the (n) value in the inter- 
val < n < 1. Phases w = — ~ (phase 1), 
?7 = | — n (phase 2), rj = 5 — § (phase 3), n = | (phase 4) exist between 
phase transition points (which is determined according to the Maxwell 
rule from the range of S'-like behaviour (between the spinodal points Zi)) 
and outside of them. At the change of the model parameter values the 



Figure 3. Field dependen- 
cy of the order parameter 
(T= 0). 



G 



regions, where metastable phases exist, can overlap, some phase transi- 
tions disappear (some intermediate phases can not be realized). In case 
1 < n < 2 the dependence r){h) is generally similar. The phase 3 and 
phase 2' at n = § — n, which now appears instead of phase 2, may play 
the role of the intermediate phases. 

On the phase diagram n — h Fig. 4 the 
thick solid line indicates the phase coex- 
istence curve and hence the possibility of 
the first order phase transition with the 
change of the longitudinal field h and/or 
electron concentration (n) takes place. 

More detail analyse of a free energy be- 
haviour shows that the above presented (in S 
,,. ,n ,. , b igurc 4. n — h phase tran- 

tms paragraph) results arc not realized. ° 

™ . ,. ,. , , , . sition diagram (T = 0). 

I he investigation of the equilibrium con- ° v ' 

ditions shows that the first order phase transition transforms into the 
phase separation. One can see regions where state with homogenous dis- 
tribution of particles is unstable (d/j,/dn < 0), and the phase separation 
into the regions with different concentrations exists (Fig. 5). The phase 
diagram n — h (Fig. 5) illustrates the separation phases. 




12 g 2g 

Figure 5. Dependence of the chemical potential n on the electron con- 
centration and n — h phase separation diagram (T = 0). 



The phase 3 splits into phase 4 (with concentration n = 0, order 
parameter n = |) and phase 1 (n = 2, r\ = — \) with weight coefficients 
1 — n and n respectively (thin dotted lines in Fig. 4 and Fig. 5). 

Therefore, more convenient and thcrmodynamically stable is the pha- 
se separated state, which is the mixture of states with different electron 
concentrations and different values of order parameters. 



1.2. Antiferroelectric type interaction. 



In the case of antiferroelectric type interaction it is convenient to in- 
troduce two kinds of sites (A-sites, B-sites). These corresponds to the 
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doubling of the lattice period [ |HJ . 

Within the framework of the MFA we shall write: 



S A S B = - V aVb + VaS z b + mS z A , (6) 

where r] A = (S%), n B = (S%). 

Then, we obtain the following expression for the model Hamiltonian: 



H 



= ~K n iAl + TliAl) + UriiAiriiAl + g(ntAl + UiAl)SiA 
-(h+ J\T\ B + J2Va)S z a , 



H 



JB 



H 



(7) 



The Hilbert space forms as a direct product of the eigenfunctions (^) for 
H A and H B operators (§) {\n. lA1 N lAl , S Z A )} © {\n lB ]N iBll Sf B )}. The 
analytical consideration in this case in general is very similar to the previ- 
ous (ferromagnetic interaction) one, but formulae are more complicated 
(see in details [ |lC|)- 

Grand canonical potential can be written in the form: 



ft = JiVAVB + y^a+ Vb) +Tln 



-Tin (l + e -K-W+ g ) + 2e -/3(-M+f)j 
-Tin (l + e -K-W-g) + 2e -' 3 (-f-i)) 



(8) 



1.2.1. [a = const regime. 




Figure 6. 
area. 



J1-J2 



The set of equations for rj A , n B is defined by the 
minimum of the grand canonical potential (^). The 
expression for the electron mean values is determined 
from the thermodynamical relation (Q). 

The form of the grand canonical potential field 
dependence Q(h) (and therefore the type and number 
of phase transitions) depends on the relation between 
parameters J\ and J2 values (Fig. 6) . There are no any specific behaviour 
when J 1 and J 2 are placed in the regions 5 and 6. The case when J\ 
and J2 are placed in the domains 1, 7, 8 is similar to ferroelectric type 
interaction. 

The location of J\ and J2 parameters within the area 4 leads to the 
possibility of the two sequential second order phase transitions: from 
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-0.8 0.0 0.8 -0.8 0.0 0.8 

Figure 7. Field dependencies of grand canonical potential within the 
2 and 4 regions respectively (T = 0). 



.1/2 
-J/2 



the ferroelectric phase (FP) with the one pseudospin mean value to the 
antiferroelectric phase (AP) and then to the FP with the another order 
parameter value (Fig. 7). 

The case when J\ and J2 belong to 
the region 2, 3 is shown in Fig. 7 and 
the corresponding phase diagram fi—h in 
Fig. 8. 

One can see that two first order phase 
transitions with the change of the field h 
and/or chemical potential a take place. 
With respect to ferroelectric type inter- 
action between pseudospins (Fig. ||) the 
phase coexistence curve is split and one 
obtains the range (the range width is equal Ji) where the AP exists 
(Fig. 8). With the temperature increase the first order phase transitions 
transform into the second order phase transitions. Hence the possibility 
of the first order phase transition from FP into AP and then the sec- 
ond order phase transition from AP into FP exist with the temperature 
increase for the narrow range of h values. 







|3> 


P> 






|4> 




ID 








/|3> 






|2> 








|T> 


|4> 




H 


-g/2 U-g/2 g/2 


U+g/2 



Figure 8. 
gram (T 



/i — h phase dia- 
= 0). 



1.2.2. n = const regime. 

As it was mentioned above, in the case of the fixed value of the elec- 
tron concentration (regime n =const) the first order phase transition 
transform into the phase separation. 

In Fig. 9 one present the phase diagram n-h when J\, Ji are placed 
in the region 2. Within the area surrounded by the lines the phase sep- 
aration into the regions with different concentrations and phases FP 
(solid lines) and AP (dotted lines) takes place. Outside of these bound- 
aries (which surround the phase separated states) the state with the 
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space homogeneity of electron concentration (FP) is stable. Between the 
boundaries one have the AP. 




Figure 9. n — h phase diagram (T = and T ^ 0). 

We would like to remind that in the /i =const regime with the tem- 
perature increase the first order phase transitions transform into the 
second order (Ji, J2 G domain 2,3). On the other hand, in the n =const 
regime this correspond to the narrowing of the range of the phase sepa- 
rated states and transform into the second order phase transition curves. 
Then the phase transition curves approach one to another and, finally, 
disappear at the certain value of temperature. 

The location of J\ , J2 within the area 4 leads to the possibility of the 
two second order phase transitions with the change of the field (similar 
to the /i =const regime). 

2. Interaction between pseudo-spins via conducting el- 
ectron. 

In the U = and Jy = limit operator ([!]) can be transformed to 
the Hamiltonian of the electron subsystem of binary alloy in the case 
of equilibrium disorder. Model ([[]) is close to the Falicov-Kimball (FK) 
model but differ in thermodynamic equilibrium conditions, i.e. in a way 
how self-consistency is achieved ({S z ) = const for the FK model and 
h = const for the PE one). 

In the present part of work we propose (for the case of the U = and 
Jij = limit) the self-consistent scheme for calculation of mean values of 
pseudospin and particle number operators, grand canonical potential as 
well as correlation functions. The approach is based on the GRPA with 
the inclusion of the mean field corrections. The possibilities of the phase 
separation and chess-board phase appearance are investigated [ . 

The calculation is performed in the strong coupling case (g t) using 
of single-site states as the basic one. The formalism of electron creation 
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(annihilation) operators a ia = b ia P^ , a ia = b ia P[~ (Pf 1 = i ± 5f ) 
acting at a site with the certain pseudospin orientation is introduced. 
Expansion of the calculated quantities in terms of electron transfer leads 
to the infinite series of terms containing the averages of the T-products 
of the Oio-, ai a operators. The evaluation of such averages is made using 
the corresponding Wick's theorem. The averages of the products of the 
projection operators P^ are expanded in semi-invariants [ 
Nonperturbated electron Green function is equal to 
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P? 



Pr 



(9) 



where e — —fi + g/2, e = —fi — g/2 are single-site energies. Single- 
electron Green function (calculated in Hubbard-I type approximation) 
is = Gh{ui n ) = (5 _1 (cj„) — ifc) and its poles determine the 

electron spectrum 



1 



1 



e 1 , u (tk) = ^(2E + t k ) ± - ^g 2 + At k {S*)g + t\ 



(10) 



In the adopted approximation the diagrammatic series for the pseu- 
dospin mean value can be presented in the form 



(s z ) 




(11) 



Here we use the following diagrammatic notations: Q — S z 
— gi(u> n ), wavy line is the Fourier transform of hopping t k . Semi-invari- 
ants are represented by ovals and contain the <5-symbols on site indexes. 
In the spirit of the traditional mean field approach [ [llj] the renormal- 
ization of the basic semi-invariant by the insertion of independent loop 
fragments is taken into account in (11). 

The analytical expression for the loop is the following: 




2 {Pti^n 

N ^ k 

n , k 

(3{ ai P+ +a 2 Pr). 



(12) 



It should be noted that within the self-consistent scheme of the 
GRPA, the chain fragments form the single-electron Green function in 
the Hubbard-I approximation and in the sequences of loop diagrams in 
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the expressions for grand canonical potential f2 and pair correlation func- 
tions ((SfSj), (Sfrij), (niUj)) the connections between any two loops by 
more than one semi-invariant are omitted. This procedure includes the 
renormalization of the higher order semi- invariants, which is similar to 
the one given by expression (pi]). 

From (|ll|) and (|l^) follows the equation for pseudospin mean value 

(5f-> = \ tanh + a 2 ({S*)) - a 1 ((^))) + In • (13) 

The grand canonical potential in the considered approximation has the 
form: 



> v - (cosh|e I (i fc ))(coshf£ II (t fe )) 

> In = 5 Y\b ) (ceo — ct\) 

N(3^ (cosh fe) (cosh §e) 



1, . f/3,. v , l+e-' 36 ! 1, . f/3, , l+e"^ 
< — (n+a2— ai)+ In ~ > +— in cosh < — 



ln COSh \ 2 (/l+tt2 ~ ai } + hl T+c^ J + ^ ln C ° Sh \ 2 k+ ln 1+e-/ 3 - 

The solution of equation (|l^) and calculation of potential fl were per- 
formed numerically. 

In the fi = const regime (in the uniform case) there exists a possibility 
of the first order phase transition with the jump of the pseudospin mean 
value and reconstruction of the electron spectrum [ |ll| . 

In the n = const regime one can see (Fig. 10) that the regions with 
dfx/dn < 0, where states with a homogeneous distribution of particles are 
unstable, exist. This corresponds to the phase separation into the states 
with different electron concentrations and pseudospin mean values. In 
the phase separated region the free energy as a function of n deflects 
up (Fig. 10) and concentrations of the separated phases are determined 
by the tangent line touch points (these points are also the points of 
binodal lines which are determined according to the Maxwell rule from 
the function /u(n)). 

The analysis of the (S z S z ) q correlator temperature behaviour shows 
that the high temperature phase become unstable with respect to fluctu- 
ations with g ^ for certain values of model parameters. The maximal 
temperature of instability is achieved for q = (tt, tt) (in the case of square 
lattice with nearest-neighbor hopping) and indicates the possibility of 
phase transition into a modulated (chess-board) phase. 

The analytical consideration of the chess-board phase within the 
framework of the presented above approximation can be performed in a 
similar way. 
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Figure 10. Dependence of the chemical potential [i on the electron con- 
centration n and deviation of the free energy from linear dependence 
AF = F(n) - [%F(2) + (1 - f )F(0)] for different T values. 



0,4 



0,2 



0,0 




From the comparison of the grand 
canonical potential values for uniform 
and chess-board phases, the (//, ft.) phase 
diagram is obtained (Fig. 11). One can 
see that chess-board phase exists as an 
intermediate one between the uniform 
phases with the different (S z ) and n 
values. The transition between differ- 
ent uniform phases (bistability) is of the 
first order (Fig. 11, dashed line), while 
the transition between the uniform and 
modulated ones is of the first (dotted 
line) or second (solid line) order. 

Appearance of various phases in the considered model remind the sit- 
uation known for the FK model with a rich phase diagram [ |l^| . However, 
contrary to this model, an existence of the phase transitions between uni- 
form phases is possible in our case. This results from the another regime 
of thermodynamic averaging (fixation of h field which is an analogous to 
the chemical potential in the FK model). 



-0,6 -0,5 -0,4 -0,3 

Figure 11. fi-h phase diag- 
ram. I - uniform phase, II - 
chess-board phase. 
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